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Abstract: The relation between finite isokinetic ther- 
mostats and infinite Hamiltonian thermostats is stud- 
ied and their equivalence in the thermodynamic limit is 
heuristically discussed. 

Studies on non equilibrium statistical mechan- 
ics progressed considerably after the introduction 
of artificial forces supposed to simulate the in- 
teraction of a "test system" with "heat reser- 
voirs", also called "thermostats". Simulations 
could be developed eliminating the need of very 
large systems to model the action of heat reser- 
voirs. The drawback is that the equations of mo- 
tion are no longer Hamiltonian. The simulations 
led to developments and to many new insights 
into nonequilibrium, particularly with regard to 
the theory of large fluctuations (fluctuation theo- 
rems, work relations and attempted applications 
to problems ranging from biophysics to fluid tur- 
bulence). An ongoing question has been, there- 
fore, whether such thermostat models are just de- 
vices to generate simulations that may have little 
to do with physical reality and, therefore, in the 
end not really relevant for Physics. There are, 
however, conjectures of equivalence between var- 
ious kinds of thermostats and the often preferred 
"infinite thermostats" which, being Hamiltonian, 
are considered more fundamental (in spite of be- 
ing infinite in size) or the stochastic thermostats. 
Here we try to substantiate, via a heuristic analy- 
sis, the equivalence conjecture between "Hamilto- 
nian" and "isokinetic" thermostats by discussing 
it in precise terms. Isokinetic will mean that ar- 
tificial forces are introduced whose role is to turn 
into an exact constant of motion the total kinetic 
energy of the particles identified as particles of 
any of the thermostats interacting with the test 
system. The novelty here is that a careful distinc- 
tion is made between the test system particles 
and the particles of the thermostats in contact 
with it but physically located in containers out- 
side the sytem (as in most real thermal baths). 
The artificial forces only act on the latter: this is 
a substantial difference from most cases consid- 
ered in the literature in which the artificial forces 
act also on the test system particles (technically 
called "bulk thermostats"): it is convenient to 
call the thermostats considered here "peripheral 
thermostats". The test system will be kept fixed 
but the thermostats will be allowed to be of arbi- 



trary size, and their behavior as the size becomes 
infinite is what will interest us. The conclusion 
is that, under a suitable assumption, a peripheral 
isokinetic thermostat becomes in the thermody- 
namic limit, when its container becomes infinite, 
completely equivalent to a Hamiltonian infinite 
thermostat: in the sense that the time evolution 
of the configurations (ie of the phase space point 
representing test and interaction systems) is, with 
probability 1, the same as that obtained by letting 
the isokinetic containers become infinite. In bulk 
thermostats there cannot be such strict equiva- 
lence because motion remains non Hamiltonian 
even in the limit of infinite systems. The analy- 
sis reinforces, as a byproduct, the identification 
(modulo an additive total time derivative) be- 
tween phase space contraction and entropy pro- 
duction. 



I. THERMOSTATS 

A classical model for nonequilibrium, for instance in 
[1], is a test system in a container fig, for instance a 
sphere of radius Rq centered at the origin O, and several 
Interaction systems containing the thermostats: we de- 
note their containers Qj and they can be thought (to 
fix ideas) as the sets flj consisting of disjoint sectors 
= G KMei > Ro,C- k, < ICk,}, J = l,.._.,n, 
kj distinct unit vectors, realized, for instance, as disjoint 
sectors in M.^, see Fig.l, i.e. as cones in with vertex 
at the origin deprived of the points inside the sphere con- 
taining the test system: for precision of language we shall 
call such containers "sperically truncated cones" ; but the 
actual shape could be rather arbitrarily changed, as it 
will appear. The terms "test" and "interaction" systems 
were introduced in [1]. The contact between test sys- 
tem and thermostats occurs only through the common 
boundaries (located on the boundary of the ball flo) of 
the test system. No scaling, of time or space, will be 
considered here. 

In the quoted reference, as well as in later related 
works, [2-4], the particles contained in r2oi ■ ■ ■ i^n were 
quantum particles and the interaction systems were in- 
finitely extended (and obeying a linear Schrodinger equa- 
tion) and each was initially in a Gibbs state at re- 
spective temperatures Ti,...,T„. Here the particles 
will be classical, with unit mass, elastically confined in 
r^o, f^i n Ar, . . . , r2n n A^ with A^ a finite ball, centered 
at O, of radius r > Rq. The temperatures in the inter- 
action systems, here called thermostats, will be defined 
by the total kinetic energies in each of them: which will 
be kept a constant of motion by adding a phenomenolog- 
ical "thermostatting force". Hence the qualification of 
isokinetic that will be given to such thermostats. More 
appropriately one should call such thermostats "periph- 
erally isokinetic" because most often in the literature the 
term isokinetic, instead, refers to systems in which the 
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total kinetic energy of all particles (in the test and inter- 
action systems) is maintained constant. The latter are 
called bulk thermostats: our models will correspond to a 
system in which no internal microscopic friction occurs 
and which exchanges energy with external systems kept 
at constant temperature. The properties that we discuss 
cannot hold for bulk thermostats. However we shall call 
our thermostats simply isokinetic except in the last sec- 
tion. The arrangement is illustrated in Fig.l below. 

Remark: peripherally isokinetic thermostats have been 
considered in the literature in simulations, [5] , and their 
physically correct behavior was immediately remarked 
sparking investigations about the equivalence problem. 
See also [6, 7]. Recently a case of a model in which only 
hard core interactions between particles were present, 
and the test system was thermostatted peripherally, has 
been studied in [8] showing the thermostat action being 
efficient and measurable even in such extreme situation. 

Phase space: Phase space Ti is the collection of locally 
finite particle configurations x = {. . . ,qi,qi, . . 

X = (Xo,Xo,Xi,Xi, . . . ,X„,X„) = (X,X) (1.1) 

with Xj c flj, hence X c f2 = U^^qQ^-, and Qi S M^; and 
in every ball B{r, O), of radius r and center at the origin 
O, fall a finite number of points oflC. 

The space 7i(Ar) will be the space of the finite con- 
figurations with X C Aj.. It will be convenient to 
imagine a configuration x as consisting of a configura- 
tion (XojXq) G H{B{Ro,0)) and by n configurations 
(X,-, X,) G n{nj n R^/B{Ro, O))), j = l,...,n. 

Interaction: The interparticle interaction will be a 
pair potential with finite range r^ and superstable in the 
sense that ip is non negative, decreasing in its range (i.e. 
"repulsive"), smooth and positive at the origin. 

Remark: Singularities like hard core could be also con- 
sidered (at the heuristic level of this paper) but are left 
out for brevity. For more general cases, like Lennard- 
Jones interparticle potentials or for modeling by external 
potentials the containers walls, see [9]. 

The potential and kinetic energies of the configura- 
tion X G H{Ar) are U{x) = Eg',g"exnA, '^(9' " 

■2 

-^(^) — SgiSXnA^ 2 ■^h^i'® the * means that the sum is 
restricted to the pairs q', q" which are either in the same 
flj or consist of two elements q' , q" of which one is in 
Qo- this means that particles in Qq interact with all the 
others but the particles in Clj interact only with the ones 
in rij U flo. The (^'s will be, for simplicity, the same for 
all pairs. 

The system in Oq interacts with the thermostats but 
the thermostats interact only with the system, see Fig.l. 




Fig.l: The 1 + n boxes Qj Ci Ar, j = 0, . . . ,n, are marked 
Co,Ci, . . . ,C„ and contain No,Ni, . . . ,Nn particles, mass 
m = 1, with positions and velocities denoted Xo, Xi, . . . , X„, 
and Xo,Xi,..., X„, respectively. The E are external, posi- 
tional, non conservative, forces; the multipliers aj are so de- 
fined that the kinetic energies Kj = ix| are exact constants 
of motion. 

Hence, if a; G W(Ar), the energy U{x) can be written 

as 



U{x) = C/o(Xo) + ^ (C/,(X^) + C/o,,(Xo, X,)) (1.2) 

and the kinetic energies will be Kjpi.j) = |X|. The 
equations of motion will be (see Fig.l) 

Xoi = -aiC/o(Xo) - ^aiC/o,,(Xo,X,) + Ei(Xo) 

Xji = -diUjCXj) - 9jl/"oj(Xo,Xj) - ajXji (1.3) 

where the first label, j = 0, . . . ,n, denotes the thermostat 
(or system) and the second the derivatives with respect 
to the coordinates of the points in the correspoding ther- 
mostat (hence the labels i in the subscipts (j, i) have 3Nj 
values); the multipliers aj are, for j = 1, . . . , n, 

^ , ^ , de/ _ ^ . _ .(Xo, X,), (1.4) 

and the "walls" {i.e. the boundaries d^i,dAr) delimit- 
ing the different containers will be supposed elastic. A 
more general model to which the analysis that follows 
also applies is in [10]. 

It is also possible to imagine thermostats acting in 
the bulk of the test system by adding a further force 
— aoXo: this is, for instance, of interest in electric con- 
duction models, [11], where the dissipation is due to en- 
ergy exchanges with oscillations ("phonons") of an un- 
derlying lattice of obstacles. Such bulk thermostatted 
systems will not be discussed because, for physical rea- 
sons, their dynamics cannot be expected to be equivalent 
to the Hamiltonian one in the strong sense that will be 
considered here. The thermostat forces would introduce 
an effective friction on the system motion not disappear- 
ing as the size of the systems grows, as it is always the 
case in bulk thermostatted systems. 
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Other thermostats considered in the Uterature, [12, 13], 
could be studied and be subject to a similar analysis, 
which would be interesting, e.g. the Nose-Hoover or the 
isoenergetic thermostats. Note that even the isoenergetic 
thermostat does not conserve Gibbs states (in presence 
of a test system). 

The equations of motion will be called isokinetically 
thermostatted because the multipliers aj are so defined 
to keep the Kj exactly constant for j > 0. The forces 
Ej(Xo) are positional nonconservative, smooth, forces. 
The numbers Nj of particles in the intial data may be 
random but will be picked with a distribution giving them 
average values of i^;'^^^ | within positive and asymptoti- 
cally Af-independent bounds as r — > oo. 



Initial data: The probability distribution /xq for the ran- 
dom choice of initial dati 
the limit as h.^ ^ oo of 



dom choice of initial data will be, if dx 11^=0 '^^^ '^^^ 



Nil 



l^o,Ar{dx) = const e ^"^^^ dx 



(1.5) 



with Ha{x) = ELo/?j(-^j(^j) - >^jNj + Uj{x)) and 



(3j , j > and /3o > arbitrary. 

Here A = (Aq, Ai, . . . A„) and T = (Tq, Ti, . . . , T„) 

arc fixed chemical potentials and temperatures {ks be- 
ing Boltzmann's constant). 

The limit /Uq as — > oo of the distribution in Eq.(1.5) 

makes sense (with particles allowed to be located in the 
infinite containers , j > 0) provided it is interpreted as 
a Gibbs distribution fiQ obtained by taking the "tcrmo- 
dynamic limit" Aj. — » oo, supposing for simplicity that 
the parameters A^ ,Tj, j > do not correspond to phase 
transition points (which would require care to consider 
boundary conditions which generate pure phases, [14]). 

It will be convenient to think always the initial data 
chosen with respect to the latter distribution: if A^ < 
DO the particles positions and velocities outside A^ will, 
however, be imagined fixed in time ("frozen", see [15, 
16]). Therefore, defining 



-l3{Kj+Uj-\N,) 



(1.6) 

and I3pj{l3,\) = limA,__oo |^^^^^| log Zj{(3, A), the ther- 
mostats density and average potential energy will be 



Un 



-knTjSj-X^Sj (1.7) 



and Sj,Uj,^kBTj will be suppposed to be the average 
density, average potential energy density and average 
kinetic energy per particle in the initial configurations: 
without loss of generality because this holds with ^q- 
probability 1 (by the no-phase-transitions assumption). 



II. DYNAMICS 

In general time evolution with the thermostatted dy- 
namics changes the measure of a volume element in phase 
space by an amount related to (but different from) the 
variation of the Liouville volume. 

Minus the change per unit time of a volume element 
measured via Eq.(1.5) is, in the sectors of phase space 
containing Nj > particles inside ArCi^lj , j = 0, 1, . . . , n, 
with kinetic energy Kj^j^^{x), 



a{x)=Y^ 



kBTj{x) 



il-{3Nj)-')+Po{Ko + Uo) (2.1) 



where kBTj{x) = ^—j^, and kBTj{x) 
oo, at least for the initial data, with fiQ probability 1. 

Remarks: (1) The dynamics given by the equations of 
motion Eq.(1.5) or by the same equations with aj = 
arc of course different. We want to study their difference. 

(2) The choice of the initial data with the distribution 
fio regarded as obtained by a thermodynamic limit of 
Eq.(1.5) rather than (more naturally) with fj.'Q j^^{dx) 



PJ' for A, 



''^'^^^'''^ = cor^st e--o(^) n - '-^) (2.2) 



dx 



with A^o, A^i, ■ 



, Nri fixed, 



dj, j > 0, and no 



lOj-nA^I "J' 

particles outside A^ is done to refer, in the following, to 
[15, 16]. A heuristic analysis would be possible also with 

this, and others, alternative choice. 

(3) The Eq.(2.2) is natural, although less convenient 
notationally, because in the case n = 1,E = and 

/3o = /3i = /3 with = fci3ri(l - 3^)^^ it is exactly 
stationary (a minor extension of [12]), if multiplied by the 

density p{x) = e ^^^j>o^^'^°''^'^ ^ which is the "missing" 
Boltzmann factor in Eq.(1.5), and therefore can be called 
an equilibrium distribution. 

Choosing initial data with the distribution /io let x — »■ 

x'^^'''°'\t)'^=Sf^'''°'-'x, a = 0, 1 be the solution of the equa- 
tions of motion with aj = {a = 0, "Hamiltonian ther- 
mostats") or aj given by Eq.(1.3) (a = 1, "isokinetic 
thermostats" ) and ignoring the particles initially outside 
Ar, [16]; and let Sj:^^x be the dynamics limA^^oo S^^^'^^^x. 

Existence of a solution to the equations of motion is 
a problem only if we wish to study the A^ ^ 00 limit, 
i.e. in the case in which the thermostats are infinite 
(thermodynamic limit). 

It is a very difficult problem even in the case in which 
a = and the evolution is Hamiltonian. For n = 1, 

def 

ai = 0, P = (3() = (3 and E = 0, a case that will be 
called equilibrium, it was shown, [15], that a solution to 
the (Hamiltonian) equations of motion exists for almost 
all initial data x chosen with a distribution obtained by 
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multiplying noidx) by an arbitrary density function p{x); 
and it is defined as the limit as ^ oo of the finitely 
many particles evolutions S^^'''^^x in f2 n A^. 

Recently, the related problem of a single infinite sys- 
tem and no thermostat forces has been solved in [15, 16] 
where it has been shown that, for a set of initial data 
which have probability 1 with respect to all distributions 
like Eq.(1.5), the Hamiltonian equations make sense and 
admit a unique solution, but the general nonequilibrium 
cases remain open. 

Therefore in the following I shall suppose, heuristically, 

a property (called below "locality of evolution") of the 
equations of motion Eq.(1.3) with and without the ther- 
mostatting forces ctj'X.j. 

The question will then be: are the two kinds of ther- 
mostats equivalent? 

This is often raised because the isokinetically thermostat- 
ted dynamics is considered "unphysical" on grounds that 
are viewed, by some, sufficient to ban isokinetic ther- 
mostats from use in physically meaningful problems, like 
their use to compute transport coefficients, [12]. The 
following heuristic considerations show that the latter 
would be too hasty a conclusion. 

III. HEURISTIC DISCUSSION AND 
EQUIVALENCE ISOKINETIC VERSUS 
HAMILTONIAN 



particle with the containers walls, afterwards they hold 
until the next collision with a new initial condition given 
by the elastic collision rule; they hold for the three dy- 
namics considered in Eq.(3.1) provided aj = in the 
second and third case and A^. is finite in the first and 
second cases. 

The first difficulty with infinite dynamics is to show 
that the speeds and the number of particles in a finite 
region of diameter r > Rq remain finite and bounded in 
terms of the region diameter (and the initial data) for all 
times or, at least, for any prefixed time interval. 

Therefore we shall suppose that the configurations 
evolve in time keeping the "same general statistical prop- 
erties" that certainly occur with probability 1 with re- 
spect to the equilibrium distributions or the distributions 
like jiQ in Eq.(1.5): i.e. density and velocity that grow 
at most logarithmically with the size of the region in 
which they are observed, [15, 16] and average kinetic en- 
ergy, average potential energy, average density having, 
asyptotically as A^. — > oo, values ^ksTj, 6j,Uj depending 
only on the thermostats parameters {Xj,Tj, j > 0), see 
Eq.(1.7). 

More precisely let the local energy in Q (1 B{^, R), ^ £ 
R^,R> Ro + r^ be 

W{x;^,R)= Yl (f + ^E^(5'-«J-)+^^r-^) 

giexnB(c,K) jj^i 

(3.3) 

with F = max \dqcp\, and its "logarithmic scale" average 



The first paper dealing with equivalence issues is [6]: W{x-^ R) 

its ideas arc taken up here, somewhat modified, and ex- ^(^) = sup sup (3-4) 

tended. A detailed comparison with [6] is in the last ?,l?l>'-v= R>r^ log ^ 

section. 

In the Hamiltonian approach the thermostats are in- and call Ho the configurations in H with 
finite systems with no thermostatting forces {aj = 0) 

the initial data are still chosen with the distribution uq „, . , 

J. , , T , tia;) < 00 and 

discussed above. Let ^ ' 



x^-\t) - (xf-^(t),xr'^(t)),.o,...,„ = sr-'>x, 



(A.,1)^ 



X^->"(i) = (Xf'-^(t),X^^-*(t)),=0,...,n = 5, 

x\t) = (XO(t),XO(i)),=o,...,n = Sf^X. (3.1) 

Then a particle {qi,qi) located at t = in, say, the j-th 
thermostat evolves, see Eq.(1.3), as 

qi{t) = qi+ I qi{t')dt' (3.2) 
Jo 

q^it) = e" /o <'o(t')dt' ^ /■* ^^„g- p,^^,r^^^„ 
Jo 

where Fi(i) = ^ag.(C/j(Xj(t)) + C/j-o(Xo(t),X,(i))). The 
above relations hold up to the first collision of the i-th 



h™ 71 7n = 0,-, 

A^^oo|ArnOj| 



lim — ^ — — - = M,-, (3.5) 



with 8 J > 0, Uj given by Eq.(1.7), if N{j, A,.), f7(j. A,,) de- 
note the number of particles and their internal potential 
energy in flj fl A^). 

The set of configurations x ^Hq has /Uq -probability 
1, [16]. 

The discussion in this paper relies on the assumptions 
1-3 below, motivated by the partial results in [15, 16], 
as it will appear shortly. It is to be expected that the 
probability distributions /x^^'^'"'*^ obtained by the 

evolution of /io with sj;^'''"''' , (a = 0, 1), and all con- 
figurations in X G Ti.Q share the following properties. 

Local dynamics assumption With ^Q-prohahility 1 for 
for X €Ho the number of collisions u{i, t, A^., a) that the 
i-th particle of x^^^''^\t') has with the containers walls 
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for <t' <t, is bounded uniformly in Ar, a, and 

(1) there is B{x,t) > 0, continuous and non decreasing 
in \t\, such that £{x^^-'°-\t)) <B{x,t), a = 0,1. 

(2) The limits x^°'\t) = liniA^^oo x'^^'-'°-\t) exist and are 
in Ho for all t, with S{x'-°\t)) < B{x,t). 

(3) x''^\t) solve the Hamiltonian equations and the latter 
admit a unique solution in Hq. 

Remarks: (a) The limits of x^^'''"'^{t), as — > oo, 
are understood in the sense that for each i the Hmits 
(^f '"^(i)) of (^f -"Hi)) exist together 
with their first two derivatives; and 
are twice continuously diffcrcntiablc in t for each i. It 
can be shown that, in the Hamiltonian case a = 0, the 
uniform bounds in (2) imply the existence of the limits, 
however they do not imply that x''^\t) G Tio, i-e. they 
do not imply the second of Eq.(3.5). 

(b) The number of points of ^''^'■''^•'(f), a = 0,1, in a 
ball B{R,i) is bounded by B{x,t)R^, for all R,S, with 
R > r^log^ and \t'\ < t. 

(c) The speed of a particle located in g' e is bounded 
by B{x,t){2log^f for \t'\<t. 

(d) Comments (b.c) say that locally the particles keep 
a finite density and reasonable energies and momentum 
distributions. 

(e) An implication is that Eq.(3.2) has a meaning with 
probability 1 on the choice of the initial data x. It is very 
important that the assumption that dynamics develops 
within Hq implies that at all times Eq.(3.5) will hold with 
Sj , Uj time independent: physically reflecting the infinite 
sizes of the thermostats whose density and enegry cannot 
change in any finite time. 

(f) The analysis of the nonequilibrium cases can be par- 
tially performed in similar Hamiltonian cases as done in 
the detailed and constructive analysis in Rcf. 16, but 
dropping the requirement in Eq.(3.5). 

(g) It seems reasonable that by the method in [16] the 
restriction of satisfying Eq.(3.5) can be removed in the 
Hamiltonian model. New ideas seem needed to obtain the 
local dynamics property in the case of the thermostatted 
dynamics. 

The multipliers a, are sums of two terms. The first is 



Applying Schwartz' inequality > exists with: 



(3.6) 



see Eq.(1.4) and the short range of the potential implies 
that the force — 9j?7oj (Xq, Xj ) is a sum of contributions 

bounded by = max |(9(/9((7)| times the number of pairs 
of particles in the band of width r^^ around the boundary 
of the container fig (because, by Eq.(3.5), £{x) < +00: 
this is of order 0{{Rlr^FS)^) if 5 is an upper bound on 
the densities near dQo- Note that such a bound exists and 
is time independent, by the local evolution hypothesis 
(above), but of course it is not uniform in the choice of 
the initial data x. 



|Xj • djUojiXo, Xj)| ^ ^ my .^ 



X| 



_J§Tj^F5_ 



(3.7) 



for Ar large and 5' = minj>o 5j, having used the first of 
Eq.(3.5). 

The second term in aj, with (ij = U{j,Ar H^j), con- 
tributes to the integrals in the exponentials Eq.(3.2) as 



2Kj 



dt" 



Uj{t)-Uj{t') 



(3.^ 



where Uj(t) is the specific energy at time t and the ~ 
reflects the use of the second equation in Eq.(3.6) to es- 
timate ^ as ksT.NUA^nn,) ' means equality up to 
quantities tending to as r 00. 

By the above hypothesis the r.h.s tends to as A^. — >■ 00 
because the configurations (initial and after evolution) 
are in Ho, hence have the same specific potential energies 
Uj (by Eq.(3.5), see also comment (e) above)), while the 
contribution to the argument of the same exponentials 
from Eq.(3.6) also tends to 1 by Eq.(3.7). 

Taking the limit of Eq.(3.2) at fixed i, this means that, 
for initial data in Ho, hence with ^0 -probability 1, the 
limit motion as Aj. ^ 00 (with f3j,Xj, j > 0, constant) 
satisfies Hamilton's equations 



qi{t) =qi+ Ut')dt', =Q^+ Fi{t")dt" (3.9) 
Jo Jo 

and the solution to such equations is unique with proba- 
bility 1. 

The conclusion is that in the thermodynamic limit the 

thermostatted evolution becomes identical, in any prefixed 
time interval, to the Hamiltonian evolution on a set of 
configurations which ha,ve probability 1 with respect to the 
initial distribution ^0, in spite of the non stationarity of 
the latter. 

In other words. Suppose that the initial data are sam- 
pled with the Gibbs distributions of the thermostats par- 
ticles (with given chemical potentials and temperatures) 
and with an arbitrary distribution for the finite system 
in 0,0 (with density with respect to the Liouville volume, 
for instance with a Gibbs distribution at temperature Tq 
and chemical potential Aq, as in Eq.(1.5)). Then, in the 
therm,odynamic limit Aj. — > oc, the time evolution is the 
sam,e that would be obtained, in the same limit, via a 
isokinetic thermostat acting in each container flj fl A^ 
to keep the total kinetic energy constant and equal to 
iNjknTj. 



IV. ENTROPY PRODUCTION 

It is important to stress that while, in the thermody- 
namic limit, the dynamics becomes the same for isoki- 



5 



6 



netic and Hamiltonian thermostats, because the thermo- 
stat force on each particle tends to 0, the phase space con- 
traction in the isokinetic dynamics does not go to zero, by 
Eq.(3.7),(3.8). Instead it becomes, up to an additive time 
derivative, see Eq.(2.1), a = Y^j^Q l^T- "^^^^ possible 
because ct is a sum of many quantities (the aj 's) each of 
which tends to in the thermodynamic hmit while their 
sum docs not. 

The interest of the remark is that Ylij>o "^jr is the 
natural definition of entropy production in both cases: 
but in the literature it is often stated (correctly so in 
the contexts) that entropy production is the phase space 
contraction, raising eyebrows because the latter vanishes 
in Hamiltonian models. 

However in finite thermostat models the phase space 
contraction rate depends on the metric used to measure 
volume in phase space: and it has been stressed that the 
ambiguity affects the phase space contraction only by 
an additive quantity which is a time derivative of some 
function on phase space. Such ambiguity will not affect 
the fluctuations of the long time averages of the phase 
space contraction which, therefore, has an intrinsic phys- 
ical meaning for this purpose, [17]. 

In both the isokinetic and Hamiltonian cases the above 
(J (which is, physically, the physical entropy production) 
differs, by a time derivative Htot, from 

-_V^/_Qj_ E(Xo) ■ Xq 

The time derivative in question here is the derivative 
of the total energy Hfot = Po{J2j>oi^3i^) + ^ji^)) + 
YjyoUoj{x)), [17]. And a generates the matter and 
heat currents, [17]. 

For this reason the equivalence conjectures, of which 
the isokinetic-Hamiltonian is a prominent example, see 
[6], [18, Sec.8],[19],[ll, Sec.6],[20],[14, Sec.9.11], to quote 
a few, are relevant for the theory of transport and es- 
tablish a connection between the fluctuation dissipation 
theorem and the fluctuation theorem, [17, 21]. 

The works [15, 16] bring the present analysis closer 
to a mathematical proof for repulsive interaction and I 
hope to show in a future work that they actually lead to 
a full proof of the locality of the dynamics, at least in 
dimension d = 1,2, for other thermostat models. 



V. COMPARISON WITH [6] & COMMENTS 

(1) Equivalence between different thermostats is widely 
studied in the literature and it is surprising that there 
are so many questions still raised about the very founda- 
tions, while little attention is devoted at trying to expand 
the analysis of the early works. A clear understanding 
of the problem was already set up in comparing isoki- 
netic, isoenergetic and Nose-Hoover bulk thermostats in 



[6], where a history of the earlier results is presented as 

well, see also [22]. 

(2) Finite thermostats acting on the boundary were stud- 
ied already in [5], in special cases, and were recognized 
to be equivalent to thermostats acting on the bulk of the 
test system. More recently, [23], isokinetic versus isoen- 
ergetic thermostats equivalence has been analyzed and 
the splitting of the phase space contraction into an en- 
tropy part and an "irrelevant" additive time derivative 
has been first stressed (see also the later [10, 24]) and 
related to the interpretation and prediction of numerical 
simulations. 

(3) The basic idea in [6] for the equivalence is that the 
multipliers defining the forces that remove the heat in 
finite thermostat models have equal average ( "equal dis- 
sipation") in the thermodynamic limit, [6, Eq.(15)]: thus 
making all evolutions equivalent. In [6] the expectation 
of observables in two thermostatted evolutions is rep- 
resented via Dyson's expansion of the respective Liou- 
ville operators starting from an equilibrium distribution: 
equivalence follows order by order in the expansion (in 
the joint thermodynamic limit and infinite time limit) 
if a mixing property, [6, Eq.(23)], of the evolution with 
respect to both the equilibrium and the stationary dis- 
tributions is assumed. The method is particularly suit- 
able for bulk thermostatted systems close to equilibrium 
where application of Dyson's expansion can be justified, 
at least in some cases, [25]. 

(4) The main difference between the present work and [6] 
is that here, even far out of equilibrium, we discuss equiv- 
alence between the boundary thermostatted dynamics 
and Hamiltonian dynamics: therefore we compare a situ- 
ation in which the average value of the dissipation (ana- 
logue of [6, Eq.(32)]) is 7^ with one in which it is 
exactly, at least formally. 

This is achieved by showing that the multipliers in the 
models in Fig.l vanish in the thermodynamic limit not 
only in average but also pointwise with probability 1; this 
is in agreement with the results in [5] and provides more 
theoretical grounds to explain them. 

It also means that in boundary thermostatted systems 
the analogue of [6, Eq.(32)] does not tend to when 
N ^ 00 although the analogue of the average of the 
multipliers, corresponding to [6, Eq.(33)], does. 

(5) In bulk thermostatted systems there cannot be equiva- 
lence between the Hamiltonian and the isokinetic dynam- 
ics in the sense discussed in this paper, i.e. identity of 
the dynamics of individual particles. However, as dis- 
cussed already in [6] , the expectation values of extensive 
observables could hold. On the other hand the analy- 
sis in [6] should be extendible to cover also the bound- 
ary thermostatted systems because, while the dissipation 
{i.e. entropy production) does not vanish in the thermo- 
dynamic limit, the average of the multipliers still does, 
sec (3) above, and this is what is really needed in [6]. 

(6) Neither Dyson's expansion convergence questions nor 
time-mixing properties, on which [6] is based, enter into 
the present analysis: but the assumptions needed on the 
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dynamics (local dynamics) are still strong and are only 

under partial control via the theory in [15, 16]. 

(7) An important question is whether taking the time 
t —> oo limit after the thermodynamic limit A^. ^ oo 
(when, therefore, the dynamics arc identical) the proba- 
bility distribution sj.^^ /iq tends to a limit /z, and /U still 
attributes probability 1 to Ho- this is an apparently 
much harder question related to the difference between 
the transient results and the, deeper, steady state results, 
[26]. 

(8) Finally: the choice, made here, of dimension 3 for 
the ambient space is not necessary for the analysis. Di- 
mension d = 1,2,3 would be equally suited. However it 
is only if the thermostats containers dimension is d = 3 
that the system with infinite thermostats is expected to 
reach a stationary state: if d = 1,2 the equalization of the 



temperatures is expected to spread from the system to 
the reservoirs and proceed indefinitely tending to estab- 
lish a constant temperature over larger and larger regions 
of size growing with a power of time, [27]. 
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